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Abstract 



I present a study of the influence of disorder in thin magnetic films on the switch- 
ing behavior of small magnetic cells of well defined shape and size. The disorder 
considered arises from randomly oriented crystalline grains of different shape, size, 
and crystalline orientation which gives rise to locally fiuctuating intrinsic anisotropy 
directions and strengths. The study comprises a theoretical investigation of a dis- 
ordered Stoner Wohlfarth model, as well as micromagnetic simulations. I show that 
the fluctuations in the total anisotropy and therefore in the switching fields are 
controlled by a single dimensionless parameter. The theoretical findings are well 
confirmed by micromagnetic simulations of many different samples. 
PACS numbers: 85.75.B, 85.70.K, 75.75 
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I. INTRODUCTION 



With the development of integrated magneto-resistive memory devices the need to pro- 
duce very many magnetic cells with reproducible magnetic switching behavior has arisen. 
As the dimensions of these devices reach orders of a few lOOnm, the switching behavior is 
to good approximation that of a single domain particle, and mainly controlled by the total 
magnetic anisotropy. The magnetic anisotropy contains contributions from various sources, 
but for such small devices the contributions from the shape anisotropy and from the crys- 
talline (i.e. material) anisotropy dominate. For a given magnetic material, the former can 
be controlled to a good degree by the aspect ratio and the thickness of the magnetic cells, 
while the latter depends on the actual microscopic structure. Typically magnetic films are 
polycrystalline with grain sizes of the order of a few nm to a few lOnm. The grains are 
oriented randomly. Since the preferred directions of the crystalline anisotropy are defined 
by the crystal axes, this leads to locally varying anisotropy contributions. Depending on the 
relative strength of shape anisotropy and crystalline anisotropy, one may therefore expect 
a more or less pronounced random component in the switching field. It is the purpose of 
this paper to assess quantitatively how much fluctuation one may expect depending on the 
various parameters in the problem. 

A lot of work in material science has been done over many years to examine how magnetic 
parameters are influenced by the composition and micro structure of thin films^^^. Models 
of disordered ferromagnets were studied analytically by Ignatchenko and co-workers in the 
late 1970s and early 1980s (see^° and references therein) as well as by Kronmiiller and co- 
workers^^. They derived "laws of approach of the magnetization to saturation". Such laws 
were used as early as 1931 for the characterization of materials^^. The laws derived by 
Ignatchenko et al. allowed to extract quantitative information about the correlations of lo- 
cal magnetic anisotropy^^ . The focus of these works was, however, on average macroscopic 
characteristics, like coercive fields or remanent magnetizations, and not on fluctuation of 
switching fields in small magnetic elements. Ref.^^ deals with the latter problem by mi- 
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cromagnetic simulations, but without a simple model presented the physical insight and 
predictive power of this analysis is limited. Also in the context of recording media, granular 
materials were extensively simulated^^. Experimentally it is very difficult to distinguish the 
contributions of different origins to fluctuations in switching fields. In particular, fiuctuations 
in the shapes of magnetic cells may easily mask fiuctuations arising from the material. 

In the following I will introduce a disordered Stoner Wohlfarth model, a simple model 
that allows to address the question of disorder induced fiuctuations of the switching fields 
analytically. I will show that within this model the joint probability distribution P{K, 7) 
for the total anisotropy strength K and the overall preferred direction 7 is controlled by a 
single dimensionless parameter u. If we denote by Ki the uniaxial shape anisotropy, by Kcr 
the crystalline (i.e. bulk single crystal) anisotropy, by N the typical number of crystallites in 
the sample, and by a a number of order unity that depends on the shape of the crystallites, 
the parameter v is given by 



The fiuctuations of the switching field decay like 1/ for large u. All of this will be derived 
in Sec. II. Sec. Ill is devoted to numerical verification by micromagnetic modeling, and Sec. IV 
contains the conclusions. 



Suppose the magnetic cell consists of crystallites, with respective volumes Vi, easy axis 
direction {I — 1, . . . , N) and crystalline anisotropy K^r. For the moment let us consider 
the simplest case in which each crystallite gives rise to a uniaxial anisotropy. Further assume 
that the magnetization M is uniform across the sample, with components uniquely specified 
by the angle 9, Mj. = M^cos^, My — M^sin^, where Mg is the saturation magnetization. 
This approximation works well for sub-//m sized magnetic cells, which are too small to hold 
domain walls. In an external magnetic field in the cell plane, H = {H^, Hy), the total energy 
density then reads 




(1) 



II. THE DISORDERED STONER WOHLFARTH MODEL 
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N 

E/V = Ki sin^ 9 + '"i - 7;) - H^M, cos 9 - HyM, sin 9 , (2) 

1=1 

where vi = Vi/V are the volume fractions of the crystaUites {V denotes the total volume 
of the magnetic cell). Note that the forms of the crystallites do not enter at this point. 
The first and the last two terms on the right hand side correspond to the ordinary Stoner 
Wohlfarth modeP^ 

In order to figure out the total anisotropy resulting from eq.(2) we need to understand 
how the different anisotropy terms add up. Let us start by looking at just two uniaxial 
anisotropies Ki and K2 with preferred angles 71 = and 72. The energy density is then 
given by 

E{9)/V = Kism^9 + K2sm^{9 -^2) ■ (3) 

We can rewrite this as 

E{9)/V ^e + Ksm^{9--i), (4) 

where e is a constant independent of 9, and K and 7 are the new anisotropy strength and 
preferred direction, respectively. By expanding the sin'^ terms in both cqs.(3) and (4) one 
easily convinces oneself that the three parameters e, K and 7 are related to Ki, K2 and 72 
by the three equations 

+ X2COS272 = Xcos27, (5) 
K2 sin 272 = K sin 27 , (6) 
K2 sin^ 72 = -ST sin^ 7 + e . (7) 

I will use the convention that all anisotropy constants are positive < Ki,K2iK < 00) and 
all uniaxial preference angles are counted in the interval — 7r/2 < 72,7 < 7r/2. The system 
of equations is then solved uniquely for all Ki,K2, and 72 by 

1 . . . /iri + i^^cos272\ 
7 = 2^ig'^(72) arccos ( — 1 , (8) 

K ^^JKl^^K^[^^2K^K^^^^2: (9) 
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and a third equation determining e, which is, however, irrelevant for the following. Note 
that the arccos function returns a value in the interval . . . tt, which by the pref ac- 
tor 1/2 is remapped to the interval 0...7r/2. Since 71 was chosen as zero, the sec- 
ond angle determines the sign of the angle of the overall preferred direction. Fig.l 
shows how the overall angle depends on 72 for various values of the ratio r = K2/ K^. 

FIGURES 
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FIG. 1. Resulting angle from adding two uniaxial anisotropies (eq.8). For r = K2/K1 < 1, the 
maximum angle reached is smaller than 7r/4. Only if K2 dominates (r > 1), the angle 7r/2 can be reached. 
The values of r chosen are 0.1 (the curve with the smallest slope at 72 = 0), 0.5, 0.7, 1.0 (a straight line 
with slope 1/2), 1.5, 2.0, 5.0, and 10.0 (almost a straight line with slope 1.0). Only the first quadrant is 
shown, as 7 is an odd function of 72. 

For K2/K1 < 1 the maximum angle reachable for the new overall preferred direction is 
smaller than 7r/4 — the Ki term always dominates and keeps the preferred direction close 
to zero. For Ki — K2, an overall angle 7r/4 can be reached at 72 = 7r/2, but eq.(9) shows 
that then the overall strength K goes to zero: With two orthogonal preferred axes with 
anisotropy of equal strength, the total anisotropy vanishes indeed. For i^2 > -^1, the overall 
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preferred angle is dominated by the second anisotropy and therefore angles beyond 7r/4 can 
be reached. 

Concerning the total strength of the anisotropy, eq.(9) tells us that K is obtained similarly 
as if two vectors with lengths and directions K2 and directions 0, 72 were added — with 
the only difference that the angle between the two vectors has to be replaced by twice 
its value in the rule for standard vector addition. Thus, the non-linear addition of uniaxial 
anisotropy terms is, for what concerns the resulting anisotropy strength, replaced by a linear 
vector addition, where the length of the vectors is given by the strengths of the anisotropy, 
and their enclosed angle is twice the angle between the original preferred directions. 

Coming back to the disordered Stoner-Wohlfarth model, eq.(2), we now understand 
that the sum of anisotropy energies with random strength and preferred directions, i.e. the 
intrinsic anisotropy = K^r YaLi vi sin^(^ — 7;), corresponds to a random walk in the plane: 
For the vector addition the relative angles are uniformly distributed in the interval [— tt, 7r[, 
and the length of step / is the random term K^rVi. The unusual addition (8) of angles not 
withstanding, the resulting final orientation of the vector in the random walk is distributed 
uniformly over [— tt, 7r[, for all since each single angle is. The random walk leads to a 
distribution of Ki that is for 1 well approximated by a Gaussian centered at — 

and with a variance cx^ = (Kf) = NK'^^ivf). If we assume that the typical volume fraction 
is given by a typical crystallite dimension ai as {vfy^'^ — yfoL{a^)ll? — ^/a/N, where a is a 
numerical prefactor depending on the distribution of crystallite areas (and thus on the shape 
of the crystallites) we obtain the scaling behavior a = ^/aKcr/ VN, and the joint probability 
distribution P{Ki,ji) of the intrinsic anisotropy and the preferred direction. 



for -7r/2 < 7i < tt/2 and < < 00. 

In order to determine the distribution of the total anisotropy, it is most convenient to 
start with an expression for the joint distribution of and the total preferred angle 7, 
employing once more the law for adding uniaxial anisotropics, eqs.(8) and (9). We have 




(10) 
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p(/^^7) 




(11) 



The two integrals are easily performed. We can then transform the distribution back to 

P{K,'y), express K with a dimensionless parameter k as K = kKi, and thus arrive at the 
final form 



with the dimensionless parameter u defined previously in eq.(l). As is now obvious, the 
total distribution of anisotropics is uniquely specified by this parameter u, and so are all 
statistical properties of the switching fields. 

The distribution is centered around k = 1 and 7 = 0, and in fact diverges for these values 
for all parameters u. Typical values of u for magnetic cells may vary over a large range. A 
rectangular 800x400x5 nm^ Permalloy cell leads to a shape anisotropy Ki — 44.4-lO^erg/cm^. 
Assuming a crystalline anisotropy Kcr = 1000 erg/cm^ and a typical crystallite size of 20nm, 
we have u — 1.58 • 10^. Larger crystallites (say {ajy^"^ — 50nm), a smaller aspect ratio, and 
thinner films may reduce this value. A rectangular 500x400x3 nm^ cell has v — 10.2 • 10^. 

The distribution P{K, 7) is most relevant for rotational remanence experiments on arrays 
of nominally identical magnetic cells^^. In these experiments a strong magnetic field is 
applied at an angle /3 relative to the nominal easy axis of the cells (as defined by shape 
anisotropy). Then the magnetic field is switched off, with the direction of the field preserved 
until zero field is reached, and one measures the remanent magnetization along the nominal 
easy axis as function of the angle p. In the strong magnetic fields (i.e. field values outside 
all astroids of the cells), a cell will always align to good approximation to the field. But 
when the field is switched off, cells that saw a positive field component along their actual 
easy axis will remain in a state magnetized along their positive actual easy axis, while those 
which saw a negative field component relative to their actual easy axis will fall into a state 
magnetized along their negative actual easy axis. Thus, all that matters is the distribution 




(12) 
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of preferred angles 7, but not the strength of the anisotropy. Let us call "x-axis"the nominal 
easy direction, the average total magnetic moment of a single cell along the x-axis at 
zero field, and n the number of cells in the array. If all cells had 7 = there would be a 
sharp jump of the total of the array at /5 = 7r/2 from = nMc to = —nMc, if (3 is 
cranked up from zero to vr. A finite width of the distribution of 7 is reflected directly in the 
width of the transition. Integrating out K in (12) we find the distribution of the preferred 
angles alone. Fig. 2 shows the result obtained from numerically integrating 

/■OO 

P(7)= / dKP{K,-f). (13) 
Jo 



Note the rather non-Gaussian profiles, in particular the pronounced cusps at zero angle, 
even for very small values of v where the distribution is almost homogeneous over the entire 
angle interval — 7r/2 . . . 7r/2. 




FIG. 2. Left: Distribution ^(7) of the preferred angle, after integrating out the anisotropy strength K. 
The parameter v ranges from v = 0.01 (no symbols, flat curve) to = 10 (triangles). Right: Scaling of the 



width of the distribution ^(7) with the parameter v. 

The standard deviation of the distribution scales as Xj^fv^ as is shown in Fig.2. For 
comparison with actual experiments it is useful to note, however, that the cells in an array 
are typically not identical, even concerning their geometry. Fluctuations in the shape due 
to lithography errors and etch processes will lead to additional fluctuations in Kx that can 
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mask fluctuations due to the intrinsic material anisotropy. 



A. Switching field fiuctuations 

Let us now have a look at the consequences for the distribution of switching fields. As 
is well known, the energy density (2) without the fluctuating term {Kcr — 0) leads to a 
stabihty region in the {H^, Hy) plane given by the astroid^^ 

A..o(9)^^!^||^ = -cos»e, (14) 

^.W-^5^ = «in'«. (15) 

The two equations are obtained by setting simultaneously dgE{6)/V = and dgE{6)/V = 0, 
which leads to the field combinations Hy{6)) where the nature of equilibrium points 

changes from stable to unstable (or from unstable to stable). A magnetic cell in a stable 
state with Mx < switches to a state > if a field combination is applied that lies 
outside of the boundary (14), (15) in the positive half plane hx > (and correspondingly 
for < 0). 

Including the random crystalline contributions in eq.(2) we now have 

E{e) /V^K siii^ie - 7) - MsHx cos 6 - M^Hy sin 6 (16) 

for the total energy density for a given cell (up to the irrelevant global shift e) , where K and 
7 are distributed according to (12). Since the easy axis direction defines the orientation of 
the coordinate frame relative to which the Stoner-Wohlfarth astroid is measured, it is clear 
even without calculation that the energy density (16) leads for each K and 7 again to an 
ideal astroid that is rotated by the angle 7 and changed in size by a factor K/ Ki (if magnetic 
fields are still measured in units of 2Ki/Mg, the normalization used in eqs.(14,15)). In the 
exphcit form of the stabihty boundary derived again from dgE{6) = = dlE{9), 

hM = ^^^^^ = (cos(3^ - 27) + 3cos(^ - 27)) , (17) 



hy{e) = ^^^^^ = (- - + 3sin(^ - 27)) , (18) 
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the rotation of the astroid is somewhat obscured by the fact that 9 is not the polar angle 
of the astroid, but just a parameter in a parametric representation. One easily convinces 
oneself, though, that acting on {hx,Oi hy^) with a rotation matrix corresponding to an angle 
7 and with an overall factor K/Ki reproduces (17), (18). 

Eqs.(17) and (18) can be used together with (12) to calculate an average astroid as well 
as fluctuations around it. It turns out, however, that for this purpose it is more convenient 
to go back one step in the calculation and keep the total intrinsic anisotropy Ki separate 
from the deterministic anisotropy Ki, since Ki is simply Gaussian distributed, eq.(lO). This 
will prove useful in Sec. II B as well, where we will look at the combination of random cubic 
crystalline anisotropy with uniaxial shape anisotropy. The latter situation will not even 
allow for defining an overall anisotropy with a single anisotropy constant. Let us therefore 
start from an energy density 

E{9) /V = Ki sin^ 9 + sin^{9 - 7^) - M,//^ cos 9 - MsHy sin 9 (19) 

with Xj, 7i distributed according to (10). If we express in units of Ki, — kiKi, we 
arrive at 

h,{9) = MfEM = -ti (cos(3e - 27,) + 3cos(^ - 27,)) - cos^ 9 , (20) 

ZKi 4 

hy{0) = ^f/^^ = (sin(3^ - 27O - 3 sin(^ - 27,)) + sin^ 9 . (21) 

Before presenting the results for the mean values and fluctuations of the switching fields, let 
me note that depending on the measurement performed different ways of averaging might 
be relevant. If we are interested in the fluctuation of the switching fleld for a given direction 
of the applied field and if the fiuctuations arc small, i.e. z/ 3> 1, the fiuctuations for a fixed 
parameter 9 are relevant. The parameter 9 will be given to good approximation for all 
astroids by the unperturbed astroid (14), (15), 

/H \^/^ 

9 = arctan . (22) 

However, when the fiuctuations are larger, or if we are interested in swept astroids where Hy 
is kept fixed and is swept, so that by definition all fluctuations are in H^, the situation 

10 



is much more complex. One will then have to calculate the relevant 9 for each realization of 
the disorder separately. This will be discussed elsewhere^'^ . 

Here I will assume that the disorder induced fluctuations are small and that the point 
on the astroid where we are interested in the variation of the switching fields is sufficiently 
well described by 9 as obtained from eq.(22). Since ki enters only hnearly in (20), (21) and 
is distributed Gaussian, and since 7^ is distributed uniformly over the interval — 7r/2 . . . 7r/2, 
we find immediately that the average astroid is the ideal astroid given in eqs.(14),(15). The 
standard deviation of the switching field in units of He — 2Ki/Ms is given by 

H 1 

crHx,Hy = -p— fV5±3cos 261, (23) 
4 y i/ 

where anx comes with the positive sign under the square root, any with the negative. Thus, 
the fluctuations in the switching field scale like l/y^- For soft magnetic materials like 
permalloy and magnetic cells with aspect ratios not too close to unity, Ki is dominated by 
shape anisotropy, which in turn is proportional to the aspect ratio over a wide range. We 
conclude that the "array quality factor" He/ (Thx for fixed cell width should be proportional 
to the aspect ratio of the cells to the power 3/2. 



B. Cubic Crystalline Anisotropy 

A material with cubic anisotropy has an energy density that depends on the direction 
cosines a^i OLy and ct^ of the magnetization with the crystal axes according to 

E/V = C{alal + alal + ay,) , (24) 

where C is the lowest order cubic anisotropy constant^^. I will assume in the following that 
the crystallites all have a z-axis perpendicular to the film, i.e. the films are supposed to 
be well ordered in 2;-direction. This is a natural assumption for flat crystallites with lateral 
dimensions of a few lOnm and films only a few nm thick, even though uniform distributions 
of the crystal axis on cones have been observed in 50nm thick films^. Furthermore, if the 
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magnetization is effectively restricted to the plane of the film (for sufficiently thin films this 
is always the case), we have = 0. I parameterize the magnetization in the plane again 
by an angle 9 with respect to the x-axis. We then have ax = cos 6*, ay = sin 6*, and the 
expression reduces to 



Cubic anisotropy projected to a (001) crystal plane thus leads to a four-fold (or bi-axial) 
symmetry in that plane. And instead of the four-fold jagged astroid for purely uniaxial 
materials, the stability curve is now eight-fold jagged. Note that the shape anisotropy of 
the cells is still uniaxial, though. The disordered Stoner Wohlfarth model for materials with 
the projected cubic crystalline anisotropy thus reads 



where I assume again that the orientations 7^ of the crystallites are distributed uniformly 
over the entire relevant interval, i.e. — 7r/4 < 7^ < 7r/4. For adding only cubic anisotropics 
the same rules (8), (9) apply as were derived from adding uniaxial anisotropics. Indeed, in 
the derivation we can just replace 9 — > 29, and replace Ki, K2 by two corresponding cubic 
anisotropy constants Ci and C2, and everything else goes through as before. The same is 
true for the random walk picture. Thus, many cubic anisotropy terms added up lead again 
to a cubic anisotropy with a distribution of the overall C, and 7 given by 2Pj(Ci,7), see 
eq. (10). The prefactor two is due to the fact that 7 now covers only half the previous 
interval. Things are different, however, when we combine the total cubic anisotropy with 
the uniaxial shape anisotropy. Obviously, the result will neither be a pure cubic anisotropy 
nor a pure uniaxial anisotropy, but rather a sum of two such terms. The boundary of the 
stabihty region derived from doE{9) = = dlE{9), with 




(25) 



E/V = Ki sin^ 9 + -Y,Vi sin^{2{9 - 7^)) - H^M^ cos 9 - HyM^ sin 9 , 



(26) 



i=l 




(27) 



now reads 
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hM = -MfEM = I (5 cos(3^ - 47,) + 3 cos(5^ - 47.)) - cos^ 6 , (28) 
ZKi o 

hy{0) = ^^^^M. = 1(5 sin(3^ - 47,) - 3 sin(5^ - 47^)) + sin^ 6 . (29) 

Depending on the relative strength Cj = Ci/Ki and orientation of the total cu- 
bic intrinsic anisotropy, this boundary may be rather different from the ideal Stoner- 
Wohlfarth astroid, as for example depicted in Fig.3. Depending on the parameters, 
little twists arise that might not always be resolvable in experiments, and give the 
impression of astroids broadened in one direction, or of kinks in the astroid sides. 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
FIG. 3. Stability boundary for a combination of cubic (crystalline) anisotropy with uniaxial (shape) 
anisotropy. Moving from the curve with largest width in y— direction to the one with smallest width in 
y direction, the curves correspond to the parameters 1.) q = 0.5, 7, = 0.5, 2.) q = 0.1, 7^ = 0.5, 3.) 
Cj = 0.1, 7i = 0, and 4.) Cj = 0.5, 7^ = 0.0, respectively. 

As is evident from eqs.(28) and (29), the average astroid for fixed parameter 6 is again 
the ideal astroid. The standard deviations of the switching fields are given by 



O'HxMy 



He 1 ; 

Vl7± 15cos2^, 



(30) 



with the plus sign pertaining to (Thx, and the minus sign to any Note that the total 



uncertainty of the switching field Jcth^ + crjjy is independent of the parameter 9 for both 
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cubic and uniaxial crystalline anisotropy. 



III. NUMERICAL SIMULATIONS 

I performed micromagnetic simulations for many different cell types, varying cell size, 
aspect ratio, thickness, and material properties. A commercially available simulation pack- 
age was used that allows to mimic a poly-crystalline structure with an adjustable average 
size of the crystallites. Each crystallite is then assigned a random preferred direction. Typ- 
ically, 10 to 100 different disorder realizations were used for each type of cell simulated, 
and average value and standard deviation of the switching fields in easy direction (actu- 
ally under a small angle of 4 degrees with respect to the easy axis in order to avoid the 
numerical problems related to catastrophic switching) were calculated. Fig.4 shows a cu- 
mulative plot for all samples with uniaxial or cubic crystalline anisotropy of the standard 
deviation of the switching field (in units of Ki) as function of the parameter For the 
calculation of v, the numerically determined average value of Ki was used, related to the 
switching field by He — 2Ki/Ms. For both types of crystalline anisotropy the decay of the 
fluctuations like 1 / is well observed over almost five orders of magnitude. For very large 
u the fluctuations seem to decay slightly slower, but they might be limited by the finite 
field resolution, as well as the intrinsic fiuctuations of the simulation program. Ideally one 
would expect all curves to collapse on a single one. The simulations show that the numeri- 
cal constant a in the definition of u does depend somewhat on the nominal sample properties. 



Note that a negative anisotropy constant corresponds to a preferred axis rotated by 90 degrees. 
As the crystallite axes are uniformly distributed over the full tt interval, P{K, 7) should not depend 
on the sign of Kcr- This was checked numerically by using some samples with negative Kcr- 
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FIG. 4. Standard deviation of the switching field fluctuations in units of the average switching field 
for various systems. Empty circles: Ellipses 200 x 400 x 4nm^ with Kcr = 1-0 • lO^erg/cm^; empty 
squares: ellipses 300 x 600 x 4nm'^, Kcr = 1-0 • lO'^erg/cm^; empty diamonds: ellipses 300 x 450 x 4nm^, 
Kcr = 1.0- 10^ erg/cm^; empty triangles up: ellipses 200 x 400 x 4nm^, Kcr = — 15- lO^erg/cm^; full triangles 
left: rectangles 200 X 600 X 4nm'^ , uniaxial ii'cr = — 15-lO^erg/cm'^; full triangles down: 200x400x4nm'^, uni- 
axial Kcr = 1.0 -10^ erg/cm^; full triangles right: ellipses 200 x 400 x 4nm^, uniaxial Kcr = — 15- lO^erg/cm^; 
pluses: ellipses 300 x 600 x 4nm^, uniaxial Kcr = 1-0 • lO^erg/cm^. The lines are guides to the eye for data 
of the same nominal system, for which the crystallite size was varied, typically between 5 or lOnm up to 
lOOnm. The dashed straight line indicates the l/v^ behavior. 



I have presented a study of the influence of crystalline disorder on the switching behavior 
of small magnetic cells. Within a Stoner-Wohlfarth model with random anisotropy contri- 
butions I have derived the joint-probability distribution of the overall anisotropy strength 
and direction. The form of the distribution implies a dependence of the switching field 
fluctuations on a single parameter u, see eq.(l), in the form of a l/v^ behavior. Also, the 
width of the transition in rotational remanence experiments should scale as 1/ y/u, and a 



IV. CONCLUSIONS 
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broadening of the transition due to crystalline disorder should lead to a rather remarkable 
line shape. Micromagnetic simulations confirmed the scaling with v both for uniaxial and 
cubic crystalline anisotropy. 
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